Evaluating management regimes for European beech forests using dynamic programming

Aim of study: This contribution describes a systematic search method for identifying optimum thinning regimes for beech forests (Fagus sylvatica L.) by using a combination of optimization heuristics and a simple whole stand growth prediction model. Area of study: Data to build the model come from standard and management forest inventories as well as yield tables from the Northern and Western part of Germany and from southern and central Denmark. Material and methods: Growth projections are made from equations to project basal area and top height. The remaining stand variables are recovered from additional equations fitted from forest inventory data or acquired from other authors. Mortality is estimated through an algorithm based on the maximum density line. The optimization routine uses a two-state dynamic programming model. Thinning type is defined by the NG index, which describes the ratio of the proportion of removed trees and basal area with respect to the same proportion before thinning. Main results: Growth equations fitted from inventory data show high goodness of fit with R values larger than 0.85 and high significance levels for the parameter estimates. The mortality algorithm converges quickly providing mortality estimates within the expected range. Research highlights: The combination of a simple growth and yield model within a Dynamic Programming framework in conjunction with NG values as indicators of thinning type yield good estimates of practical thinning schedules compared to thinning recommendations provided by diverse authors.